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Abstract. Geometrically, foams or covalent graphs can be decomposed into successive layers or strata. 
Disorder of the underlying structure imposes a characteristic roughening of the layers. Our main results 
are hysteresis and convergence in the layer sequences. 

1) If the direction of construction is reversed, the layers are different in the up and down sequences 
(irreversibility); nevertheless, under suitable but non-restrictive conditions, the layers come back, exactly, 
to the initial profile, a hysteresis phenomenon. 

2) Layer sequences based on different initial conditions (e.g. different starting cells) converge, at least in the 
cylindrical geometry. Jogs in layers may be represented as pairs of opposite dislocations, moving erratically 
due to the disorder of the underlying structure and ending up annihilating when colliding. 

PACS. 83.70.Hq Heterogeneous liquids: suspensions, dispersions, emulsions, foams, etc. - 61.43.-j Disor- 
dered solids - 87.18.Bb Multicellular phenomena: computer simulations - 68.35.Ct Interface structure and 
roughness 



1 Introduction 



Despite the broad range of their different material reali- 
sations (liquid froths, metallic microstructures, polymeric 
foams, etc.), foams and analogous disordered structures 
have strong structural similarities . Even if universal- 
ity has not been firmly established for these structures, 
many features obviously do not depend on the details of 
the constituents, of the force fields, etc. nor on the precise 
values of metric quantities. Whence, our choice of describ- 
ing foams and random patterns at the level of topology. 

To account for correlations and statistics beyond the 
one body properties, it is natural to define configurations 
in terms of topological distance. Indeed, in foams, nearest 
neighbour cells are clearly defined by sharing an interface. 
In the dual network, the cells are represented by points 
connected by bonds, one for each facet in real space. Thus 
the dual is agraph, which is sufficient to define topological 
distance [|,||]|,§. 

Covalent structures, like those occurring in glassy ma- 
terials, are also described in terms of graphs: the atoms, 
molecules or clusters sit at the nodes and the covalent 
links define the bonds. Since covalent interactions are car- 
ried by quantum electron clouds, defining the bonds is not 
always free from ambiguities. This is still more so when 
non covalent interactions are involved. In these cases, ge- 
ometric constructions, such as the Dirichlct-Voronoi one, 
may complement or replace physico-chemical criteria. 



The minimal number of (dual, in the case of foams) 
bonds, or steps, needed to join two nodes defines a topo- 
logical distance (or simply distance, since no other notion 
of distance will be considered here). 

A layer is the set of nodes / cells at a fixed distance 
j from an origin O. Partitioning the entire foam into suc- 
cessive layers j = 1, 2, 3, . . . provides a stratification of the 
cellular pattern. 

There are many reasons for improving our understand- 
ing of stratifications. The number of nodes / cells, in suc- 
cessive layers — the population, for short — gives almost 
the same information as the pair correlation function 
||,[7]]. As is well known, the correlations are related to the 
response of the system to all kinds of solicitations. In dis- 
ordered materials, this question is of particular interest: 
is the response coded in the geometry and how ? Re- 
ciprocally, beyond elasticity, external actions may modify 
the structure. How ? Aging is a common characteristic of 
glassy materials which almost never reach equilibrium; ag- 
ing may occur spontaneously or under external influence, 
often in an inhomogeneous way. 

All these questions involve structure. Our purpose, here, 
is to analyse some of the fundamental geometric tools, 
to set the ground for further research. Ultimately, energy 
should be considered. But, in complex systems, the step 
from geometry to energy is often easier than understand- 
ing the geometry. Foams are paradigmatic in this respect: 
to first approximation, energy is film length (in 2D) or 
area (in 3D) times a constant (surface tension). 
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Viewed as dynamical processes, considering j as time, 
layer sequences represent the successive stages of signals, 
fronts, epidemics propagating at unit velocity. There is 
a close analogy with aggregation-deposition and related 
problems (cf. j?J and refs. therein). 

One of the differences, however, is that the underlying 
foam is given a priori in its full integrity. The stratification 
is a supplementary structure — an ordered partition — of 
the foam. It is therefore necessary to disentangle what 
is general to layers, from what depends specifically on 
the underlying structure. Notably, the same structure can 
have many possible stratifications. A preliminary answer 
to the question of sorting these stratifications will be found 
here through convergence. 

The arbitrariness comes from the choice of the origin. 
It corresponds to the large collection of different stratifi- 
cations that are possible for partitioning a given foam or 
covalent structure. What symmetry covers these equiva- 
lent choices ? For example, in 0, the leading asymptotic 
behaviour of the layer population Kj was found (numer- 
ically) to be independent of the central cell (the origin). 
Here, we give a strong support to this hypothesis by show- 
ing that the stratifications do actually converge. The cen- 
tral cell may even be replaced, as the origin, by a central 
set of cells. 

Another distinction is of importance: the elementary 
models of aggregation, such as the Eden model, are ran- 
dom processes on regular lattices || . In our case, the un- 
derlying structure — a foam — is random, with quenched 
disorder, whereas the process — stratifying — obeys a fixed, 
deterministic, rule (without any randomness). Some geo- 
metrical features, such as roughness, are similar in both 
types of systems J?J. In the present paper, we insist on 
aspects which are more specific to the second class: (irre- 
versibility, hysteresis, reciprocity, convergence. 



1.1 Layers 

In the foam, we classify the cells in terms of topological 
distance: layer number j (layer j or lay(j') for short) is the 
set of cells at distance j from O. The origin O may be a 
single cell or a cluster of cells (concentric geometry) or a 
connected set such as a row of cells (going once around 
the cylinder in cylindrical geometry, or infinitely long in 
open Euclidean geometry). 

(In general, the dual of a cellular structure is a multi- 
graph. For almost all natural foams, it is a simple graph, 
that is, neighbouring cells have at most one facet in com- 
mon. This technical assumption is not essential, and it 
could easily be lifted if necessary, but it makes the pre- 
sentation simpler.) 

For cellular structures, this defines layers through the 
dual. It is also possible to operate in the direct cellular 
network as follows |],[|. Consider an initial cluster O, 
called the origin and labelled j = 0; the cells in con- 
tact with O constitute the first layer. Then, inductively 
for j = 1,2,.. ., layer j is made of all the cells, not yet 
counted, which are in contact with layer j — 1. 



If, as in modelling chemical structures, the origin is 
a single vertex — an atom in the compound — , then the 
layers are coordination shells |l^,|ll|,|2). So layers, coronas 
[ [13| and coordination shells are synonyms. We also name 
them strata, because they partition the foam — the set of 
vertices in the dual — into an ordered collection of subsets, 
making altogether a stratification (or foliation). 

The embedding space is either Euclidean (the plane in 
2D) or a cylinder equivalent to a domain of bounded base 
with periodic boundary conditions in the x direction (a 
circle in 2D) and infinite along the axis of the cylinder (y 
coordinate) . 

In layer j, 

— every cell is neighbour of at least one cell in lay(j — 1); 

— some cells, called regular, are also neighbour of cells in 

lay(i + i); 

— the other cells, not sharing any edge (facet in 3D) with 
cells in lay(j + 1), are called defects. 

The first statement above is part of the definition of 
layer j; the next two are definitions of regular and defect 
cells in the layer. As we shall see, defects are sources of 
frustration, curvature and non-triviality of the stratifica- 
tion. 

In summary, a stratification £ = {£j}j>o is a partition 
of the foam (or of the set of nodes in covalent graphs) into 
layers — the strata £j, j — 0, 1, 2, . . .. Each layer is the 
set of cells at distance j from O: £j = {c | dist(c, O) = j}. 

Shell hj is defined as the outer boundary of layer j; it 
is the contour lying between layers j and j + 1. 



1.2 The columnar model 

This toy model (the columns) is a useful laboratory for 
the structure of foams and as a model of growth. It is a 
lattice version of the Poisson partition of Fortes |l4j . We 
use it for illustration but most of the features presented 
here are valid generally, not limited to this example. 

The relevance and limits of this model were discussed 
in Q. For completeness' sake, we recall the basics here. 

The model is a 2D packing of columnar cells, each of 
width 1 (in the horizontal, x, direction) and of random 
length s (height in the vertical, y, direction). The sizes (s 
is both length and area) of the individual cells are taken as 
independent random even numbers, identically distributed 
with exponential law : 

Pr(s) = — - z s / 2 , s = 2,4,6.... (1) 

z 

The parameter z has a fixed value in ]0, 1[. It controls the 
mean cell size through (s) = 2/(1 — z). Unless otherwise 
stated, we will take z = 1/2, (s) = 4. 

The foam lies on a semi-infinite vertical cylinder, mean- 
ing that it is periodic, with period L, in the x direction. 
The height s of each cell is an even random number. With 
ground yo(x) — x mod 2 (crenellated profile), this ensures 
that the vertices have coordination 3, as in real foams. 
The system is unbounded in the positive y direction. 



C. Oguey et al.: Stratifications of cellular patterns 



3 




To avoid overloading the pictures, the graphical con- 
vention of Figure [l] (right) will be used: the cell bound- 
aries are not drawn; the lines are layer boundaries (= shells 
hj,j = 0, 1, 2, 3, . . .). The top square of each cell is marked 
by a dot (.) when the cell is regular, by a cross (+) when 
it is a defect J|,||,0. 

2 Up and down: irreversibility 

We compare different stratifications on the same cellular 
pattern. In this section, we build two sets of layers, one 
with distance increasing upwards (stratification from the 
bottom up), the other with distance increasing downwards 
(stratification from an origin at the top). Later, in Section 
0, we compare stratifications rising in the same direction 
but based on different origins (grounds). The question is 
whether they match, and, if so, how ? 

2.1 Up and down 

Starting from an origin Ao = O or a ground ho, if we 
build the layers upwards A = {Ai, A%, ■ ■ ■ , Aj, . . .}, stop 



at some j = d and then, taking layer Ad as a new origin O' 
(equivalent to setting shell hd-i as a starting profile h' ), 
build new layers A' = {A' — O', A[, A' 2 , . . .} downwards, 
this new stratification A' does not coincide with the former 
(even if we compare just the regular parts). The top most 
layer of A coincides with the origin of A' , by construction. 
But then, some cells, qualified as defects in upward layers 
switch to being regular in downward layers and vice versa, 
etc. Since these switches cumulate during buildup of the 
stratification, we may expect that the coherence between 
the two reverse stratifications A and A' is rapidly lost. 
This appears to be the case, at first sight (see Figure || 
left and middle). 

Notably, the last shell going down, h' d _ 1 , is different 
from the upward starting ground ho. It lies in a neigh- 
bourhood of ho, but it is different. This difference will be 
described in Section 

2.2 Back up 

What happens if we go back again ? Take the last down 
shell h' d _ x as a new ground, h® , and build another strat- 
ification B = {Bo = A' d , B\, B2, ■ ■ ■} climbing up again. 
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15 20 

, 15, up. 



5 10 15 

Idem, descending. 



Fig. 1. L = 30 : stratifications of 15 layers up, down and up again. The ground j 
hj(x) of the two upwards stratifications are identical for j > 12. 



5 10 15 

Idem, up again. 
is highlighted in each case. The profiles 



This third stratification B — {Bj} is different from the 
first two: different from A' because of irreversibility; differ- 
ent from A because the new ground, Hq , is not, in general, 
a shell hj of the first stratification (Figure || right) . 

Nevertheless, after climbing up, building B up to Bd, 
the last profile fits exactly the same profile as the first 
crest: h^_-y = hd-i, Bd = Ad- This will be proved later. 

Further up and down processes repeat A' and B. In- 
deed, the next stratification B' (downwards) is degenerate 
with A' since it starts from the same origin, and so on. 
Recall that all these stratifications are based on the same, 
fixed, but random, foam. 



2.3 Hysteresis 

We have therefore a hysteresis cycle, caused by the pres- 
ence of defects. Indeed, defect-free stratifications are re- 
versible. Examples of these are the rows and columns 
allel to the square basis in le Caer's construction 
the vertical columns in the columnar model, or even the 
horizontal layers {^} after defect coalescence. 



All these are flat or pur e gauge models like the Mattis 
model for spin glasses JIt]]. But, let us stress this point, 
these models admitting defect-free stratifications are not 
generic. Notably, in le Caer's model, there are special cor- 
relations between neighbouring cells |jl6). Topologically 
random foams are not flat. Hysteresis might even be taken 
as a measure of non trivial disorder. 



3 The geometry of layers 

In this section, we present the stratifications as analogous 
to foliated structures. In particular, a proof is given of the 
fact that the extreme layers are exactly recovered by the 
down-up procedure. 



3.1 Layers as sets 

The distance between sets A, B is defined as 

dist(^, B) = min{dist(a, b)\a G A, b E B}. (2) 
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In this sense, layer j, as a set of cells, is at distance j from 
the origin O which can consist of more than a single cell. 
In fact, by definition, all the cells of lay(j) are at distance 
j from O. But the converse is not true: not all the cells 
of O are at minimal distance j from lay(j). This is a first 
indication of irreversibility. 

3.2 Geodesic sections 

Let us go on. By definition, any cell Cj of lay (j) has at least 
one neighbour Cj-i in lay(j — 1); Cj—i has a neighbour Cj_2 
in lay(j — 2) etc. down to some root cell cq in O. Thus, 
there is always at least one connected section linking any 
cell Cj in layer j to O. In the dual, these connected sections 
are lines of minimal length j = dist(cj, O), i.e. topological 
geodesies, linking O and lay(j). Moreover, these sections 
consist of regular cells exclusively. In particular, defects in 
O cannot be root cells. 

Stratification is analogous to foliation in differential 
geometry. The layers are the leaves, and any section can 
serve as base space (isomorphic Z or some subinterval of 
Z). 

The layer structure is robust along these sections: there 
is exactly one cell per layer crossed. Along these lines, 
each step is a move from a layer to the next one — up- 
ward or downward. Thus, the sequence of layer numbers 
j = 1,2,... coincides with topological distance along these 
lines (counted, respectively, from bottom up, or from top 
down ) . The set of linking geodesies constitutes an orthog- 
onal skeleton for the stratification. 

As already stated, there is a section linked to every cell 
in lay(j), but not every cell o of O is at distance j from 
lay(j); only the root cells are. Linking geodesies starting 
from different top cells may fuse on the way down. So the 
whole set (of linking geodesies) is a forest with branches 
attached to every cell of lay(j) but only a few root cells 
in lay(O) = O. (Note that there may be more than one 
geodesic connecting twoeiven cells). Another forest pat- 
tern was introduced in || . 

The space left — that is, the part of the foam not cov- 
ered by the skeleton — is the place where irreversibility 
occurs; the layers down differ from the layers up. 

3.3 Up and down revisited: parallel layers 

Two sets A and B are parallel if there is a positive number 
d such that i) all the a € A are at the same distance d to 
B and ii) all b G B are at the same distance d to A. 

With respect to stratifications, where the sets are sets 
of cells and distance is topological distance, parallel sets 
enjoy special properties. If A and B are parallel at distance 
d: 

— In the stratification based on A, B is a subset of the 
dth layer: lay(O) = A => B c lay(d). 

— All the cells of A are root cells. 

— Conversely, A is a subset of the dth layer based on B: 
lay(O) = B =s> A c lay(d). 



— All the cells of B are root cells in this stratification. 

In the notations of Section |2[ we now show that the 
layers Ad and A' d are parallel at distance d. The proof 
requires a few basic (in)equalities. 

Up? 

In the up stratification Aq = O, .., Aj, .. — as in any stratif- 
ication — dist(lay(j), O) = j only implies dist(lay(j), o) > 
j for an arbitrary cell o of O. Equality holds if and only 
if o is a root cell Co for some geodesic section. Moreover, 
equality must hold for at least one cell; there always is at 
least one root cell in O. 



Down | 

Consider the down stratification A' , A[, A' 2 , ■ ■ . , A' d . A' Q 
consists of all the cells of Ad of the up stratification. Call 
{hi} the cells of layer A' d . {bi} includes all the root cells of 

0. All the others must lye below {bi} because they satisfy 
strict inequality: dist(o, Ad) > d. 

Parallelism J 

Thus, A' and A' d are parallel. 

Indeed, i) dist(6, A' Q ) = d, Vb G A' d holds by definition 
of layer A' d . To see that ii) dist(a,v4^) = d, Va G A' , note 
that the construction of the layers implies dist(a, A' d ) > d. 
On the other hand, in the up stratification, there is a cell 
Co G O at dist(a, Co) = d; being root, eg also belongs to 
A' d . Therefore dist(a, A' d ) < d, which proves equality ii). 

Remarks 

1. Note that A' , which was set equal to Ad, contains no 
defect for the downward stratification. Indeed, any cell 
c G Ad is at distance 1 of a (regular) cell of Ad-i and 
any regular cell of A^—i is reached this way, implying 
A 1 d ^ 1 C A[ . Now a defect Cd in A' would be at distance 
at least 2 from A' ly in contradiction with Cd G Ad = 
{cldist^A^Hl}. 

2. Going up and down establishes a 'reciprocity' relation 
between layers A' = Ad and A' d , slightly stronger than 
parallelism. Such a reciprocity does not hold in general; 
most often, two layers in the same sequence are not 
even parallel. In order to get a pair of reciprocal layers, 
a precise procedure must be followed, such as the up- 
down trick. 

3. Apart from that, nothing special is assumed on either 
the foam or the original layer. The point where we turn 
back (j — d) is chosen arbitrarily; layer Ad is abso- 
lutely normal, with neither more nor less defects than 
any other. Actually, if d is large enough, the system 
basically forgets its initial conditions. Another strati- 
fication, in the same sense but starting from another 
origin, would generate the same layers at sufficiently 
large j. 
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4. A stack of successive layers has minimal thickness when 
it is delimited by a pair of reciprocal layers. With 
the previous notations, this means dist(o, A4) > d = 
dist(c, Ad) for all o in Aq, c in A' d , whatever O we start 
from. 

5. Reciprocity does not mean reversibility of the process 
(layer sequence) nor does it coincide with parallelism. 
It implies parallelism, but parallelism is weaker be- 
cause it misses some completeness condition, as shown 
in the following example. 



Convergence 

From numerical simulations on columns and topological 
foams with randomly generated h§, we observe that the 
stratifications {Aj},{Bj} converge: for any h^, there are 
integers J, k such that Bj = Aj+k for a ll j > J- 

The rate of convergence will be discussed later (Sec. 
4.3). First, we analyse the phenomenon in terms of dislo- 
cations. 



An example: concentric stratifications 

Take a single cell O = {0} as origin and build the concen- 
tric stratification around it. Then o must be a root cell. 
Going up and down (out and in, implying a new strati- 
fication inwards) brings one back to a cluster C contain- 
ing the starting cell o possibly surrounded by other cells. 
Lay(j), which, in this case, is the topological circle of ra- 
dius j, and its centre {0} are parallel (according to our 
definition). But they are not 'reciprocal'; up-down does 
not come back to only {o}. Cluster C, on the other hand, 
is both parallel to, and in reciprocity relation with, the 
topological circle since it was constructed so. 

Notice that o is parallel to any topological circle around 
it (any azimuthal layer at distance j = 1, 2, . . .). But clus- 
ter C is parallel only to some specific circle (s), where the 
turn back is done, or could be done, in order to get C 
exactly. 

Irreversibility, or hysteresis, is the fact that, in between 
C and the circle j, the outwards and inwards stratifica- 
tions are different; this is visible only if j > 2. 



4 Convergence of the stratifications, 
dependence on ground 

The choice of the origin O is arbitrary. One may choose 
a single cell, and obtain concentric layers. But choosing 
a horizontal ground is better adapted to cylindrical ge- 
ometry. Consider a definite foam on a cylinder. Call A — 
{A }j>o the stratification based on y = Hq (x) ~ (a con- 
nected set of cell boundaries; y is the coordinate along the 
cylinder axis). 

For the same foam, we could take as origin another 
profile {}Iq[x) I x = 0, ..,L — 1} following other cell edges. 
Let B = {Bj}j> be the stratification based on h$ . How 
do A and B compare ? 



Global shift 

If h,Q is a shell of A, say = h£ for some integer k, then, 
trivially, Bj = A, J+ k, Vj > 0. The two layer sequences 
are identical; only their label differ by an integer k (an 
irrelevant phase shift). Therefore, only profiles ho with 
centre of mass near y = need be considered. 



4.1 Dislocation pairs in the stratifications 

Apart from the flat ground Hq(x) — 0, the simplest start- 
ing ground is a 'podium': 1%q{x) = 1 for x+ < x < X-, 
= otherwise (in vertical units of layers). The steps at 
x+, X- are a pair of dislocations in the stratification, with 
strengths +1, -1. Because of periodic boundary conditions, 
the strengths must sum up to 0. 

Let us compare A, based on a fixed ground, with a 
stratifications B, based on a podium tig of width w and 
of height 1 in units of ^4-layer thickness. At j = the 
dislocations are at the ends of the podium: x+(0) and 
X-(0) with |x+(0) — x_(0)| = w. Choosing the maximal 
distance w ~ L/2 will give an estimate of the convergence 
time for more general situations. 

At any later 'time' j, away from the dislocations x±(j), 
the two layer systems (profiles, inclusions, etc.) are the 
same, except for a shift of 1 in numbering between the 
two dislocations. The differences arc confined to the region 
near x + and x_ where the numbering makes steps. 

As can be seen in Figure [|, where are only marked cells 
which are defects in one stratification but not in the other, 
the differences look like two random walks which ulti- 
mately annihilate, as in a 'diffusion-reaction' phenomenon. 

The convergence occurs at time of first collision J, 
when the opposite dislocation meet for the first time and 
cancel. The layers agree from there on because, for a fixed 
underlying foam, the process 

. . . -> lay(j - 1) -> lay(j) -> lay(j + 1) -> . . . 

is deterministic. 

Due to periodicity in the x direction, the dislocations 
may fuse on one side (with vanishing — > region), 
or on the other (vanishing — = region). Conver- 
gence means that Aj — Bj in the former case, Bj = Aj+\ 
in the latter, for j > J. 

Incidentally, in a crystalline foam, the analogous tra- 
jectories would be periodic in space (ballistic regime). There- 
fore convergence would occur in time j = J linear in L (or 
not at all, when the lines cc+ , x_ are parallel) . 

In random foams, convergence depends on disorder. 
Here it is faster than in standard diffusion; the spreading 
grows wi th t ime to the exponent 1/z = 2/3 instead of 1/2. 
(see Sec. [If). 

Remark The various stratifications are made over a given, 
random structure. Drawing the successive layers is there- 
fore an entirely deterministic process over the same 
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Fig. 3. Two examples of layer convergence. For a given, 
columnar foam, we compare two stratifications: B is based 
on a podium at j = 0; the reference A is based on a flat 
ground. Only the cells which are defect in one stratification 
(A or B) but not in the other are marked. The podium has 
width w = L/2, and centre at L/2 in one case (+) and at 
L/2 + 6 in the other (x). In both cases, the pair of disloca- 
tions annihilates at some time j — J (different in each case). 
In the first case, the layers end up in phase. In the second case, 
the final time shift is one (as if the podium had covered a full 
layer). The sample contains (L — 100) x 400 cells. 



\ 



Fig. 4. (Ahj) = (mint, \hf — hf +k \) as a function of j for 
L = 100, 200, 400, 800. (Average over 50 stratifications B x 50 
different foam structures). Fits are of the form cexp(— j/r). 
For each foam, the 50 stratifications B have randomly chosen 
grounds. Each stratification contains more than 10000 layers. 




Fig. 5. Correlation time r in function of L, log-log. The line 
is 0.162 x L 15 . 



random structure. Convergence is like many of these 
mechanisms for finding successive key cards in a given, 
shuffled pack. Once two stratifications are in phase at 
some time J, they remain in phase thereafter. 

Because of periodic boundary conditions, a general 
ground ho can always be decomposed into dislocation pairs 
(+/-1 steps). When, initially, there is a large density of dis- 
locations (highly corrugated )), many dislocation pairs 
cancel at small j because the partners are initially close 
to each other; this holds for random (diffusion) and crys- 
talline (ballistic) foams. The ultimate convergence of the 
stratifications is controlled by the few dislocations that 
sur vive at longer time (j). This is further analysed in Sec. 
4.3. 



4.2 Attractor 

Clearly, the outcome of the convergence is a layer system 
— a stratification — which is a stable attractor. 



Specifically, there are two stationary stratifications: 
one up and one down. 

4.3 Convergence rate 

If the random motion of the dislocations is governed by 
some cooperative phenomena related to roughening, then 
J should be the time needed to reach £j = L, £ being the 
correlation length along the layer. 

First, at a fixed sample width L, convergence occurs 
at an exponential rate. This has been checked by measur- 
ing the mean distance between the profiles, (Ahj), where 
Ahj = minfe \hf — h£ +k \ (Figure^). Indeed, the correlation 
"time" r, defined by (Ahj) oc exp(— j/r), is finite as long 
as the maximal possible distance between dislocations is 
bounded, as it is for finite L. 

The long time pseudo-diffusion process is manifest in 
the dependence t(L) on sample size L. In the columnar 
model, which has been shown to fall into the KPZ uni- 
versality class @,|l^l, the characteristic "time" r ~ (J) 
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scales as r 



-1.5. 



1.5 is the dynamic exponent. 5.2 Anisotropy 



This prediction is well confirmed by our simulations. In 
the range of large L, the plot (Figure |J) shows a scaling 
behaviour fitting a power law r ~ L in accordance with 
KPZ. 



5 Discussion, conclusions, perspectives... 

5.1 Summary 

For foams or random covalent structures, we have shown 
that the layer sequences are irreversible. The stratifica- 
tions in one direction and the other differ even if the two 
sequences share a whole layer (this can be done by the up- 
down trick). The hysteresis between up and down strati- 
fications is due to topological defects, inherently present 
when the disorder in non trivial. The up-down procedure 
leads to topologically parallel layers at distance d, enjoy- 
ing special properties. 

— Reciprocity: each one may be reached from the other 
by building a sequence of d layers. 

— Minimal thickness of the enclosed stack: any set of d 
successive layers ending at one of the parallel layers, 
but based on another initial condition at j = 0, will 
have a thickness larger than the strip bounded by the 
parallel layers. 

On a given fixed foam, the stratifications based on dif- 
ferent origins converge to an attractor, one for each of 
the two directions (in cylindrical geometry). This pair of 
attractive stratification appears to be specific of the un- 
derlying cellular pattern. The characteristic time for con- 
vergence, t ~ L 15 , agrees with KPZ universality class, 
as long as the probability distribution decays rapidly for 
cells with a large number of sides n (exponentially, or as 
n~ K with k large enough) |19|j . This has been confirmed 
by numerical simulations on the columnar model. 

As a consequence, stratifications built on two foam 
samples differing only by local perturbations — topological 
transformations like neighbour exchanges, cell birth or co- 
alescence, etc. — will also converge, eve n if t he convergence 
is, in practical respects, slow (see Sec. 4.3). 

As proven in Sec. ^, the first set of properties — hysteresis, 
reciprocity, minimal thickness — hold generally, for any 
type of foam or graph. 

Convergence and attractors, however, are still conjec- 
tural. They essentially follow from an interplay between 
determinism of the process and randomness of the land- 
scape. Simulations of rectangular foams with periodic bound- 
ary conditions in one direction, infinite in the other direc- 
tion, confirmed the phenomenon and gave us quantitative 
results on the rate of convergence, its scaling properties 
and its relation to roughness. 

The main biases of our model are that the disorder is 
confined to one direction and that the width of the sys- 
tem is finite. These two aspects, local and global, deserve 
separate discussions. 



We think that the columnar nature of our model has neg- 
ligible influence on our observations; our conclusions hold 
more generally. Convergence was probed in the disordered 
(vertical) direction, where randomness provides a good im- 
itation of more realistic foams. 

Preliminary simulations of topological foams — generated 
by operating a large number of randomly distributed topo- 
logical transformations as in [||[20| — show the same prop- 
erties as those observed in the rectangular model: hys- 
teresis, of course, but also, to some extent, convergence of 
stratifications, etc. 

Notice that convergence is observable and measurable 
in any type of foam, not only columnar. This is a signifi- 
cant improvement with respect to (tJ, where most of the 
analysis was based on height h(x), which is rather specific 
to the columnar model. 



5.3 Boundary conditions 

The extension to foams in other types of spaces is twofold. 

As already argued, and shown on an example in con- 
centric geometry, parallelism, irreversibility and hystere- 
sis, which can be tested in finite regions, occur quite gen- 
erally: in planar or 3D foams, embedded in Euclidean or 
curved spaces. 

As to convergence, it holds unambiguously only in cylin- 
drical foams. These boundary conditions introduce a def- 
inite length-scale into the system. For quantum gravity, 
this might be an unbearable hypothesis. At more common 
scales, cylindrical geometry is quite frequent. Condensed 
matter, zoology or botanic, etc, provide lots of examples 
with tubules, channels, stems, stalks, straws,... 

When dealing with other boundary conditions, the ques- 
tion of convergence is not straightforward. There are el- 
ementary obstructions to the onset of a uniform conver- 
gence in concentric geometry. However, convergence may 
still be true in a weaker sense, either in the mean over 
each layer, or restricted to sectors of prescribed aperture. 
All these questions are under current investigations. 
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